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Abstract. A self-similar solution for the 3D axi-symmetric radiative MHD equations, which revisits the formation and 
acceleration of accretion-powered jets in AGNs and microquasars, is presented. The model relies primarily on electromagnetic 
extraction of rotational energy from the disk plasma and forming a geometrically thin super-Keplerian layer between the 
disk and the overlying corona. The outflowing plasma in this layer is dissipative, two-temperature, virial-hot, advective 
and electron-proton dominated. The innermost part of the disk in this model is turbulent-free, sub-Keplerian rotating and 
advective-dominated. This part ceases to radiate as a standard disk, and most of the accretion energy is converted into magnetic 
and kinetic energies that go into powering the jet. The corresponding luminosities of these turbulent-truncated disks are 
discussed. 

In the case of a spinning black hole accreting at low accretion rates, the Blandford-Znajek process is found to modify the total 
power of the jet, depending on the accretion rate. 
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1. Introduction 
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£j , Recent theoretical and observational efforts to uncover the mechanisms underlyi ng jet formation in AGNs and micro-quasars 
^ ' leave little doubt a bout their linkage to the accretion phenomena via disks (Mir abel 200 it iLivio 1999k iBland ford 2001: 
• • l lHuieiratetal. 20031) . 

. £h ' innermost part of a disk surrounding a black hole is most likely threaded by a large scale poloidal magnetic field (-PMF) of 
, external origin JBlandford & Payne 1982). Unlike standard accretion disks in which the sub-equipartition magnetic field (-MF) 
^ ' generates turbulence that is subsequently dissipated and radiated away, a strong PMF in excess of thermal equipartition may sup- 
press the generation of turbulence and acts mainly to convert the shear energy into magnetic energy. A part of the shear-generated 
toroidal magnetic field (-TMF) undergoes magnetic reconnection in the launching layer, whereas the other part is advected out- 
wards to form relativistic outflows (Hujeir at et al. 2003I) . 

Indeed, the total power of the jet in M87 is approximately Lj m 10 44 erg jBicknill jfelSegehnaji^^. This is of the same 
order as the power resulting from a disk accreting at the rate M ~ 1.6 x 10~ 3 AlEdd <Di Matteo et a. 2003). Similarly, the vari- 
abilities associated with the microquasar GRS 1915-105, as have been classified by (Belloni et al. 20001) . consists of a state "C" 
which corresponds to the case when the innermost part of the disk diminish jFender 19991) . probably a phase in which the super- 
luminal jet is fed with rotational and magnetic energies. 

Another two additional basic questions related to jet formation are: 1) What are the contents of these jets? specifically, are they 
made of electron-positron (e~/e + ), electron-proton (e~/ p) or a mixture of both. 2) How significant is the role of the central object 
in the formation and acceleration of jet-plasmas? 

It is generally accepted that electron temperature § e ( = kT/m e c 2 ) of order unity can be easily achieved in accretion flows onto 
BHs, and specifically in the vicinity of the event horizon. At such high temperatures pair creation becomes possible. The ef- 
ficiency of pair creation depends mainly on the interplay between the optical dept h to scattering T n _ pk^ c H, on the so ca lled 
"compactness" parameter C — a?oq~H 2 and on the strength of the magnetic field (Bio rnsson & Svensson 199HlEsin 19991) . H 
here denotes the thickness of the disk or a relevant length scale, oq is a constant coefficient, p is the density and k sc is the opacity 
due Thomson scattering, and q~ is the radiative emissivity. Generally, the efficiency of pair creation decreases strongly with 



Send offprint requests to: A. Hujeirat, e-mail: hu j eirat@mpia-hd . mpg . de 



2 



A. Hujeirat and R. Blandford: formation of accretion-powered jets in radio galaxies 



Please ask A. Hujeirat: 
huj eirat @ mpia-hd . mpg . de 



Fig. 1. A schematic illustration of the jet-disk connection model. In the outer disk, r > r tr , MFs are weak and advected around by 
fluid motions. Balbus-Hawley and Parker-instabilities in combination with reconnection and fast inwards fluid-motions amplify 
the MFs and reach thermal equipartition at r tr . Interior to r tr , MFs become of large scale topology relative to the flow in the disk, 
hence adopting the profile r~ 2 to conserve the poloidal magnetic flux. Here, torsional Alfven waves are the dominant angular 
momentum carrier. They transport angular momentum from the disk upwards and form a super- Keplerian thin layer, where the 
plasma rotates super-Keplerian, virial-heated by reconnection, and start to centrifugal-accelerate outwards. 

increasing the optical depth r p . It decreases also if the strength of the magnetic field is increased. Esin (1999) has studied pair 
creation in steady hot two-temperature accretion flows for a large variety of parameters. The ions here are preferentially heated by 
turbulent dissipation and cool via Coulomb interaction with the electrons, which in turn are subject to various cooling processes 
such as Bremsstrahlung and Synchrotron emission. Esin (1999) showed that two-temperature accretion flows are almost pair- 
free, whereas single-temperature accretion flows are more appropriate for pair-dominated plasmas. However, enhancing thermal 
coupling between electrons and protons, or heating electrons and protons at an equal rate, and relaxing the stationarity condition 
may lead to e~/e + -dominated plasma, provided the accretion rate is sufficiently low. 

In the case the BH is rotating and accreting at low rates, formation of <?~/ e + jets though ele ctromagnetic extraction of rotational 
energy from the hole via Blandford-Znajek (1977) process is very likelv dRees et al. 19821) . Indeed , recent observation s reveals 
that the jet-plasma in M87 is probably an e~ / e + -dominated jet ( Revn oldsetld^9'9c)l) . Furthermore. llWardle etal. 1998T) reported 
about the detection of circularly polarized radio emission from the jet in the quasar 3C 279. The circular polarization produced 
by Faraday conversion requires the radiative energy distribution to extend down to low energy bands, which indicate that e~/e + - 
plasma population might be a significant fraction of the jet-contents. 

On the other hand, recently it has been argued that unless the poloidal magnetic field threading the event horizon is exception- 
ally strong, the total available power through th e Blandford-Znajek process is dominated by the accretion power of the the disk 
(Ghosh & Abra mowicz 19971 iLivio et al. 19991) . We note, however, that the P MF can be sufficient ly strong, if the innermost part 
of the disk collapses dynamically, while conserving poloidal magnetic flux jHuieirat et al. 20031) . and which may significantly 
enhance the efficiency of the Blandford-Znajek process. 

In this paper we present a theoretical model for the formation of accretion-powered jets in AGNs and microquasars. The model 
relies mainly on electromagnetic extraction of rotational energy from the disk plasma and deposit it into a geometrically thin 
layer between the disk and the overlying corona (Sec. 3). In Section 4 we discuss the content of the plasma in the TL. The role 
of the Blandford-Znajek process in combination with the luminosity available from the transition layer is discussed is discussed 
in Sec. 5. We compute the total luminosity of truncated disks in Sec. 6, and end up with summary and conclusions in Section 7. 

2. The model problem 

Let r = r tI be a transition radius, where the effects of magnetic fields on the dynamics of accretion flows become significant. 
This may occur if the initially weak magnetic fields in the dissipative-dominated disk-plasma are amplified via dynamo action, in 
which Balbus-Hawley (1991) instability in combination with the Parker instability results in a topological change of the magnetic 
field from a locally disordered into well-ordered large scale magnetic field (see Fig. 1). At this radius the MFs are in equipartition 
with the thermal energy of the plasma, i.e., B = E mag /E^ = 1. 

Interior to r tr the poloidal magnetic field B p is predominantly of large scale, and magnetic flux conservation implies that PMF 
increases inwards obeying the power law 1 B p ~ r~ 2 . Such a strong PMF suppresses the generation of turbulence, whenever/? 
exceeds unity. In this case, the heating mechanisms that hold the plasma against runaway cooling are magnetic reconnection, 
adiabatic compression and other non-local sources such as radiative reflection, Comptonization and conduction of heat flux from 
the surrounding hot media. 

Given a plasma threaded by an r -2 - PMF, the angular velocity must deviate significantly from its corresponding Keplerian profile. 
This is because torsional Alfven waves (-TAWs) extract rotational energy from the disk plasma on the time scale: ttaw = H^/V^, 
which decreases strongly inwards. V V A here is the Alfven speed due to the PMF. The TAW crossing time of the disk, i.e., ttaw, is 
of the same order, or it can be even shorter than the dynamical time scale Td yn = r/V^. In this case the angular velocity becomes 
sub-Keplerian, and the disk becomes pre-dominantly advection-dominated. If tjaw < T"d yn , accretion may terminate completely. 
TAWs, which carry with it rotational energy from the disk-plasma, may succeed to propagate through the overlying corona 

1 O = 2nr 2 B p ss 2nr 2 B e = const. 
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without being dissipated. This corresponds to magnetic braking of the disk, in which the rotational energy is transported from the 
disk and deposited into the far ISM. However, since the PMF-lines rotate faster as the equator is approached, the PMF winds up, 
intersect and subsequently reconnect, and so inevitably terminate the propagation of TAW into higher latitudes. Alternatively, the 
vertical propagation of the TWAs may terminate via magnetic reconnection at the surface of the disk, establishing thereby a highly 
dissipative transition layer (TL), i.e., chromosphere, where the toroidal magnetic fields intersect, reconnect and subsequently heat 
and particle-accelerate the plasma in the TL. The latter possibility is more plausible, mainly because: 

1. The vertical profile of the angular velocity generate a toroidal magnetic field of opposite signs, i.e., anti-parallel toroidal 
flux tubes, which give rise to magnetic reconnection of the toroidal magnetic field lines (-TMF). The reconnection process 
terminates the vertical propagation of the TAWs, inducing thereby a magnetic trapping of the rotational energy in the TL. As 
a consequence, a centrifugally induced outflow is initiated, which carries with it rotational, thermal and magnetic energies. 

2. Unlike stell ar corona e that are heated from below, BH-coronae have been found to be dynamically unstable to heat conduction 
( Huiei rat et al. 20 02). The effect of conduction is to transport heat from the innermost hot layers into the outer cool envelopes. 
Furthermore, noting that the Lorenz and centrifugal forces in an axi-symmetric flow attain minimum in the vicinity of the 
rotational axis, we conclude that this region is inappropriate for initiating outflows. 

3. In a stable flow-configuration, an isolated Keplerian-rotating particles in the disk region move along trajectories with mini- 
mum total energy. If these particles are forced to move to higher latitudes while conserving their angular momentum, then 
their rotation velocity becomes super-Keplerian, and therefore start to accelerate outwards as they emerge from the disk sur- 
face. Since the speed of these motions and the associated strength of the generated toroidal flux tubes increase inwards, the 
magnetic flux tubes are likely to interact and subsequently to reconnect. 

4. Observations reveals that the radio luminosity in the vicinity of the nucleus of various AGNs with jets acquire a significant 
fraction of the total luminosity. This indicates the necessity for an efficient heating mechanism that allow the electrons to 
continuously emit Synchrotron radiation on the Alfven wave crossing time. Noting that B v ~ r~ 2 , and that the plasma 
in the innermost part of the disk is turbulent-free, we propose that magnetic reconnection of the TMF is most reasonable 
mechanism for heating the plasma in the TL. This agrees with the proposal of (Ogilvie & Livio 2001) who argued that jet 
launching requires thermal assistance. 



2.1. The velocity field 

The problem we are addressing in this paper is: assume that the innermost part of the disk is threaded by r -2 - PMF. What are 
the most reasonable power-law distributions of the other variables in both, the disk and in the TL, that give rise to inflow-outflow 
configuration, and which simultaneously satisfies the set of the steady 3D axi-symmetric radiative MHD equations. As we shall 
see in the next sections, there are several signatures that hint to Q ~ r~ 5 ' 4 as the optimal profile (see Fig. 2). Nevertheless, we 
will consider this profile as an assumption. 
Therefore, in the disk region the following profiles are assumed: 

£2(r) = Q„X- 5/4 , Bp(r) = B e (r) * B X-\ (1) 
where X — r/r tr . 

To assure a smooth matching of the variables interior to r tl with those of the standard disk exterior to r tl , we set Qg = GMj r\ v 
and By = B 2 (r = r tr ) = (\2ifR„l 'fi(y - l))pT\ r=ru , Inserting £2(r) in Equation A. 2 (see Appendix), and taking into account that 
the internal and magnetic energies of the plasma at r tr are negligibly small compared to the gravitational energy of the flow, the 
radial velocity then reads: 

U 6 if) = \Ul+\-d--\)] 112 . (2) 
e r tr a 

As in classical disks, we set Uq = U(r = r tr ) = -(2/3)v/ftr, and v = a s;i HdV s for describing turbulent viscosity. 
The radial accretion rate in the disk region is obtained by 6— integrating the continuity (Eq. A.l): 

-Al" = 2rcos0 d (pV)| ed , (3) 
Or 

where yVt^ = rE^Ud is the radial accretion rate in the disk region, and Ed = r cos Qp&dO. 9d is the angle such that tan 9d = Hd/r, 
where Hd is the vertical scale height of the disk (see Fig. 3). 

Furthermore, it is reasonable to assume that the radial dependence of the density at the surface of the disk does not differ 
significantly from the density at the equator. Thus, p(r, 9 = 9d)/p(r, 9 = 0) = const. « 1. 

Similar to the disk region, the following angular velocity and the PMF profiles are set to govern the plasma in the launching 
region: 

Q(T) = Q*F~ 5/4 , B V (Y) * B r = B^Y~ 2 , (4) 
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Fig. 2. The profiles of the angular velocity in the innermost part of the disk and in the launching region. TAWs extract rotational 
energy from the disk-plasma efficiently and deposit it into the plasma adjusting to the disk-surface. As a consequence, the inflow 
adopts the sub-Keplerian profile Q ~ r~ 5 ^ 4 in the disk-region, while the outflow adopts the super-Keplerian profile Q ~ r~ 5 ^ 4 in 
the TL. Note that the up-stream boundary conditions determine whether the flow is sub- or super-Keplerian rotating. 




9 = (Equator) 

Fig. 3. The material fluxes governing an arbitrary volume cell in the TL are shown. 6w,d are latitudes that correspond to the 
locations of the interface between the disk and the TL, and to the upper surface of the TL, respectively. F z , and F R , denote 

' > r j in,out in, out 

the material fluxes across the vertical and radial surfaces. The width of the cell is Hw, and the vertical material flux across 0\y is 
set to vanish. 



where Y — r/r^h, and tbl is the innermost radius of the TL, where the effective gravity is set to vanish. This implies that at Y — 1, 
the angular velocity approaches the Keplerian values: £2* = Q W (T = 1) = ^GM/r^ L . Inspection of Eq. A. 2 implies that the 
radial velocity in the super-Keplerian region has the profile: 

t/w(r) = [t/ W o+4-^ 1 --F] 1/2 - < 5 > 

where Uw,o is the radial velocity at tbl- However, since the effective gravity has a saddle point at tbl, the radial velocity must 
vanish. Interior to Y = 1, the plasma is gravitationally bound, so that the plasma and the associated rotational energy will be 
advected inwards and disappear in the hole. In the case of a spinning BH, tbl is likely to be located inside the ergosphere. 
The outward-oriented material flux in the TL is: 

M? = rS w f/ w , (6) 
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where 2 W = J^" r cos Op^dO. 9w is the angle such that tan6*w - tan6*d = H w /r, and H w is the vertical scale height of the 
launching layer. 

Global mass conservation of the plasma both in the disk and in TL regions requires that: 
d .-,d ~ d ,-,w 

— M r = 2— M, , (7) 
Or Or 

where mass loss from both sides of the disk has been taken into account, and where have assumed a vanishing vertical flux of 
matter across # w (see Fig. 2 ). 



2.2. The geometrical thickness of the launching region 

Since the plasma in the TL rotates super-Keplerian, the horizontal component of the centrifugal force, i.e., (tan 8 / e 2 )(V 2 / V) (see 
Eq. A. 3), acts to compress the plasma in the launching region toward the equator and its collapse. Inspection of Eq. A. 3 reveals 
that there are four forces that may oppose collapse: gas and turbulent pressures, poloidal and toroidal magnetic fields. In terms 
of Eq. A. 3, the dominant terms required for establishing such an equilibria are: 



tan OVl 



^W P s^ : ion pressure (/) 

-\§§Pim ■ turbulent pressure (//) 

\\f g Bl: poloidal MFs (///) 

ll 4jBl: toroidal MFs (IV). 



(8) 



V p r de"T 

Equivalently, the relative geometrical thickness of the TL reads: 



r 



e x (Vs/V<p) : ion pressure (I) 

e x (Vtur/V^) : turbulent pressure (II) 

e x (V V JV V ) : poloidal MFs (III) { > 

e x (Vl/Vy) : toroidal, MFs (IV), 



were we have excluded magnetic tension from our consideration, because the plasma in the TL is dissipative. Global magnetic 
flux conservation, i.e., V • B = 0, requires that B p must bend as it emerges from the disk region, provided that B p is of large scale 
topology and is advected inwards with the matter. In this case, the above-mentioned third possibility, i.e., case III, applies to those 
MFs in which B x decreases vertically and therefore may hold the plasma in the TL against vertical collapse. This means that the 
MF-lines in the TL must point toward the central BH as they emerge from the disk (see Fig. 4). Although such configurations are 
not in view with our classical expectation about magnetic-induced launching, the radiative MHD calculations do not appear to 
exclude such solutions, provided that TL-plasma is highly dissipative. On the other hand, if the poloidal MF-lines in the TL point 
away as they emerge from the disk, which is the case we are considering in this study, B x must increase in the vertical direction, 
and therefore may enhance the collapse. Therefore, taking into account that the electron thermal energy and the PMF-energy are 
relatively small compared to the rotational energy in the super-Keplerian layer, and that Z? T -profile has turning points in the TL, 
we may conclude that P tur is the most reasonable force to oppose the TL-collapse. 

To elaborate this point we note that when TAWs start propagating from the disk upwards while carrying angular momentum, 
a profile such as shown in Fig. 5 would result. This profile has VQ of mixed signs, which therefore induces the formation of 
toroidal flux tubes whose the associated electric currents j move anti-parallel around the axis of rotation, and which gives rise to 
partial magnetic reconnection (Fig. 5). 

The strong shear in the TL is capable to produce magnetic energy that is in equipartition with the rotational energy, i.e., ~ V v , 
which is necessary to maintain the steadiness stationarity of the solution of the Bj - V v coupled system. 

We note that magnetic reconnection is a common phenomena in astrophysical MHD-flows, which is generally associated with 

changes in the magnetic field topology and partial loss of the magnetic flux. In the solar case, for example, magnetic reconnection 

is one of the main mechanisms underlying the eruption phenomena of the solar flares. Approximately 10 32 ergs are liberated in 

each event and on the time scale of few minutes. Although the sun is not as compact as old stars such as a white dwarf or 

a neutron star, solar flares appear to be associated with y-ray emission, electron-positron annihilation at 5 1 1 Kev, and most 

impressively, with the neutron capture by the hydrogen nucleus which occurs at about two mega electron volts, revealing thereby 

that acceleration of particles up to relativistic velocities is an important ingredient of the reconnection process. 

Yet, the TL is located deep in the gravitational well of the central BH, and if reconnection occurs, then the energetics associated 

must be much powerful than in the solar flares, and a significant populations of e~/e + are to be expected in the TL. 

In plasma physics, reconnection occurs normally on microscopic scales. In the present study we assume that reconnection can be 

described by the magnetic diffusivity: 



Vmag - #wVmr, 



(10) 
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where Z/w is the width of the TL and V tul is the velocity that is associated with magnetic reconnection. We adopt the Petscheck- 
scenario for describing the the reconnection velocity jPriest 19941) : 

Vrec = V tuI = a mag V J A /logRe, (11) 

where a mag is a constant less than unity and log(Re) is the logarithm of the magnetic Reynolds number. 
Inserting V tul = a lmg V A / log(Re) in Eq. 9/11, we obtain 2 : 

#W _ V tm - _ a mag 

~~~V^~ \og{Re) 

From the definition of the magnetic Reynolds number: 

Inertial force Bp V» Bp log Re 

Re = — - = (—)(—) = (—)( ), (13) 

Viscous force Bj V tm - Bj a mag 

we obtain an equation for the magnetic Reynolds number: 

\og(Re) - [a m , g A]Re = 0. (14) 
B p 

Obviously, a maa must be relatively small for this equation to provide a positive Re (see Fig. 6), which means that the TL must be 
geometrically thin, and that e ~ a ma „/ log Re = Hw/r. In order to assure that the outflow leaving the system has sufficient TMF 
to collimate and reach large Lorentz factors, we require that only a small fraction of the total generated TMF dissipate in the 
TL through reconnection events, while the rest is advected outwards. Therefore, the advection time scale should be equal to the 
amplification time scale of the TMF, i.e., r a d v = r amp . This yields: B V /B T = H w /r = e. 



2.3. The profile of the toroidal MF 

In the disk region, where the plasma is turbulence-free, TAWs are the dominant angular momentum carrier. Assuming the plasma 
to be instantly incompressible, the equation describing the propagation of the magnetic torsional wave reads approximately: 

^B T * V 2 A AB T . (15) 

A denotes the two-dimensional Poisson operator in spherical geometry, and Vp (= B p / yfp) is the Alfven speed. The action of 
these waves is to extract angular momentum from the disk to higher latitudes (see Fig. 4). Note that B p determines uniquely 
the speed of propagation, hence the efficiency-dependence of angular momentum transport on the Z? p -topology. These torsional 
Alfven waves (-TAWs) propagate in the vertical direction on the time scale: 

_ H d _ H A yfp^ 

t"taw — — — ■ 

Va B e 

Obviously, since attains a minimum value and Va attains a maximum value at the innermost boundary, ttaw can be extremely 
short so that, in the absence of flux losses, a complete magnetic-induced termination of accretion should not be excluded. 
To maintain dynamical stability of the disk and teh TL, extraction of angular momentum should be compensated by radial 
advection from the outer adjusting layers, i.e., 

15, dB T 
r z or o8 

where I — r 2 cos 2 pQ. is the angular momentum. 

When ^-integrating this equation from -0<j to 6>d, we obtain: 

1 c) 

--(Myn) = B g B T , (16) 
r z or 



where Bj corresponds to the values in the TL, and which can be determined by matching both the disk- and the TL-solutions. 
Taking into account that Ai r varies slowl; 
depends weakly on the radius (see Fig. 9). 



Taking into account that yVtf varies slowly with radius, and inserting D. ~ r 5 ^ 4 and B p ~ r 2 , we obtain Bj ~ r ^ 4 , which 



2 Without considering e; it will be determined later. 
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Fig. 4. Possible configurations of converging and diverging poloidal magnetic field in the disk and the TL. The left configuration 
corresponds to the case in which the magnetic pressure at the interface between the disk and the TL increases vertically, and 
therefore enhances the collapse of the TL. The right case corresponds to magnetic pressure that opposes collapse. 

2.4. The profiles of the other variables 

1. To have stationary solutions, the energy due to the toroidal magnetic field in the TL must be is in equipartition with the 
rotational energy, i.e., £^ ag = E rot . This means that: 

Vl = = V v , => P = A 2 = const. = p%. 
VP y<p 

Thus the disk must supply the plasma in the TL with the optimal material flux across the interface, so to maintain the density 
constant. 

When r-integrating the equation of mass conservation, Eq. 7, we obtain: 

■ d ■ d ■ W ■ W ■ d • W 

M r = AC tt + 2 ( M, - AC lr ) = M r =r, r + 2AM . (17) 
This gives the density profile along the equator: 

AL d AL d 
rH d U d er z Ud 

which is a function of radius and of p^. 
The density in the TL then reads: 

pw,o = (A^ - Mtrm)'*" 7/ V- 1 /2,, where T = [Vr- 1] 1/2 (19) 

Although p w o can be treated as an input parameter, we provide here its value relative to the central density at the equator. 
Taking into account that U^i vanishes at r BL and applying mass conservation to an adjusting volume cell in the TL, we obtain 
the rough estimate: 

Pw.o = [pej[&)][(^-] = [epj[e 2 ][l/e] = e 2 p d (r = r BL ). (20) 

U nw 

In deriving this result, we have used Vg d = e 2 Ud, which relies on asymptotic expansion of the variables q = qo + eqi + 
e 2 q2 + where qo is the value corresponding to an equilibrium state, i.e., hydrostatic equilibrium in the vertical direction 
(Regev & Huieirat 1988, see the references therein). Thus, replacing d M r jdr by M r /r in Eq. 7, we may obtain a more 
accurate value for the vertical velocity across the interface at 9^. 

p^r^T 

V(r,e d ) = — — . (21) 

[Mo-2pW(M w + r^)] 

The outward-oriented material flux compared to the accretion rate is displayed in Fig. 7. The corresponding 2D profiles of the 
poloidal and toroidal components of the velocity field are shown in Fig. 8. 

2. The profile of the magnetic diffusivity in the TL reads: v mag = ffwVT = erBj/ yp ~ r 3 ^ 4 . 

3. The ion temperature is obtained by requiring that the advection and heating time scales to be equal, i.e., T a d v = Th ea t, where 
Theat = C\rpT[/<^. O is the dissipation function: v mag |V x B\ 2 w v^gBj/H^. This gives T{ ~ r~ l l 2 , or more accurately, 

Ti = T°Y- l/2 , where T° = l/e 2 . (22) 
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Fig. 5. A schematic description for a time-sequence of the angular velocity profile in the vertical direction. TAWs transport 
angular momentum vertically, enforcing the matter in the disk to rotate sub-Keplerian, and super-Keplerian in the TL (top). 
The resulting Q-profiles acquire negative and positive spatial derivatives (middle), which in the presence of a large scale PMF, 
generate toroidal flux tubes of opposite signs (bottom) that subsequently reconnect and annihilate. 



The electron temperature is extremely sensitive to the strength of the magnetic field. Since the heat liberated through magnetic 
reconnection can heat both the electrons and protons equally, the temperature of the electrons can be found by equalizing 
the Synchrotron cooling rate to the heating rate. In order to calculate the Synchrotron cooling, we must first find the critical 
frequency, below which the media becomes optically thick to Synchrotron radiation, i.e, v c , below which the emission follow 
the Raleigh-Jeans blackbody emissivity profile. This requires solving the equation: 



B v dS, (23) 



where s v and B v are the Synchrotron and Raleigh-Jeans blackbody emissivities, respectively. Having calculated v c , the 
Synchrotron cooling A syn is then calculated by dividing the total luminosity from the surface of the TL divided by its cor- 
responding TOforne^vintegrated frequency from v = up to v c , where the medium is assumed to be self-absorbed for v < v c . 
This gives (Esi n etal. 19961) : 



2nkT e vl 
3H w c 2 



A S yn * TTT^r- (24) 



Synchrotron emission at higher frequencies is a ssumed to be negligibly sm all. In terms of the density and magnetic fields, 

the following approximation due to can be used (Shap iro & Teukolskv 19831) : 

A syn = 2.16 x 10 4 ^r 9 2 p 10 , where B 3 = B/10 3 G, T 9 = T/10 9 K andpio = p/lO~ w gr. 
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Therefore, from the equalizing the heating to cooling rate, we obtain the following profile for the electron temperature: 
r e = 7^' y~ 5/8 , and T° = 1/e. 

5. It should be stressed here that only a small fraction of the total generated TMF is allowed to undergo magnetic reconnection 
in the TL. This is an essential requirement for not obtaining radio luminosity that dominate the total power emerging from 
jet-bases, and to assure that the outflow is associated with TMF-energy sufficient enough for collimating the outflows into 
jets. Thus, comparing the rate of heating via magnetic reconnection with the generation rate of the TMF, we obtain: 

(1/TDiss) T gen , B T V mag V2a mag 

(1/Tgen) TDiss Bp H^V V log(fle) 

where we have used v mag = H^/V 1 ^ and B ? /B T = H w /r = e. Consequently, an e-fraction of the total generated toroidal 
magnetic energy undergoes magnetic reconnection, while the rest is advected with the relativistic outflow. 

We now turn to find out the appropriate profiles of the electron- and proton-temperatures in the disk region. The electrons in the 



Fig. 6. The profiles of aRe and log Re versus the magnetic Reynolds number Re, where a = awg^)- The intersection point Re c 
corresponds to the solution of Equation 14. 

disk region are heated mainly by adiabatic compression, conduction and by other non-local energy sources. However, they may 
suffer of an extensive cooling due Synchrotron emission. We assume that at the interface between the disk and the TL, there is 
a conductive flux that is sufficiently large so to compensate the cooling of the disk-electrons through Synchrotron emission. In 
other wards, we require that rs yn = Theat- In this case the electron temperature adopts the adiabatic profile; T e ~ r . 

How does the ion-temperature correlate with that of the electrons in the disk region? 
Assume that the elec trons and the protons to have Maxwell velocity distributions. The time scale required to establish an equi- 



where Tg = T/ IQ 9 K, and the proton number density is given as «i4 = n/10 14 . Therefore, in the disk region Ttherm is of the same 
order as or greater than the hydrodynamical time scale if n > 10 10 , which implies that T e * T p for reasonable accretion rates. In 
the transition layer, however, n(TL) ~ e 2 n(disk). Therefore T e and T p can be significantly different. 

3. Electron-positron versus electron-proton jets 

Since the flow in the TL is heated through magnetic reconnection, it is naturally to ask whether it is <?~/<? + - or e~/ /^-dominated 
plasma. 

In view of the solar flares, we may expect the plasma in the TL to accommodate different plasma-populations such as electron- 
positron, electron-proton in a see of high y-rays. The relativistic acceleration that particles experience during reconnection 
events in the vicinity of th e BH gi ve rise to different types of pair crea tion, in particular y - y, y - e T , e T - e T , y - p, e T - p 
JSvensson 1982tlBiornsson & Svensson 199lHWhite &Lightman 1 989). Among other effects, pair fraction depends crucially on 
two important parameters: the optical depth to scattering and on the so called compactness parameter. The latter parameter is a 
fundamental scaling quantity for relativistic plasmas having § = kT e /m e c 2 is of order unity. 

Taking into account that the plasma in the innermost region of the disk is freely falling, and that Ma = 2H&r < p& > U, we 
obtain that the optical depth to scattering in the innermost part of disk is: 




librium ( S pitzer 19561) is: 

Ttherm = 4.1 l(r 9 3/2 /«14)(20/Z«A) S, 



(26) 




Ks Md 

2rU 



= i 8 . 86( J_)-1/2(_£L), 



7"LSO Afedd 



(27) 
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Fig. 7. The radial profiles of the accretion rate in the disk region, and the outwards-oriented material flux in the launching layer, 
^outflow, which applies for r < r tI . Note that ^outflow vanishes at the inner boundary tbl, where the effective gravity is set to vanish. 
Assuming e = 0.1 = H^/r, the maximum material flux associated with the outflow is approximately 5% of the total accretion 
rate through the outer boundary. Thus, F out fl ow depends on the disk temperature; colder disks yield lower ,F out fl ow and vice versa. 

where tlso is the last stable radius. Noting that pw = £ 2 pd, and that the optical depth obeys a similar relation, i.e., = e 2 r|?, we 
may conclude that the TL is optically thin to scattering for most accretion rate typical for AGNs. Further, since the flow in the 
TL is steady and two-temperature, our models fall in the left-down corner of the r s - M diagram of in Fig. 10 of Esin (1999), 
from which we conclude that jets are formed in ion/proton-dominated plasma. 




4. Extraction of rotational energies from the BH and from the disk 

A spinning black hole is surrounded by the so called ergosphere, inside which no static observer is possible; any frame of ref- 
erence must be dragged by the spin of the black hole. If the plasma in the ergosphere is threaded by an external and ordered 
magnetic field, then the MF-lines must be dragged as well, thereby generating an electric field E — V x B. In the ergosphere, 
however, negative energy orbits are possible (i.e., the total energy including the rest mass of the particle is negative). If a PMF 
threads the event horizon, it can put particles on negative energy orbits. These particles will be swallowed by the BH, while 
other particles with positive energy will emerge that carry electromagnetic energy flux from the hole. These extracted ener- 
gies from the hole is used then to accelerate the plasma-particles relativistically, possibly forming e~/e + -dominated jets (see 
Blandford & Znai ek 1977tlRees e t al. 1982; Begelm an et al. 1 984). The direction and magnitude of the electromagnetic flux are 
described by the Poynting vector: 

P = -— ExB = - — VxBxB. (28) 

Using spherical geometry, the components of the Poynting flux read: 

P = -i{ (V^Bt - V r B 2 - V r B 2 g + V e B r B e ), 

{V x B r B e - V g B 2 - V g B 2 T - V v Bj), (29) 
(VgB e B r - V^BgBj - V V B 2 + V t B t B T )}. 

In the case of a spinning BH surrounded by a plasma with zero poloidal motion, i.e., V = V r = Vg = 0, a PMF of external 
origin in the vicinity of the event horizon is predominantly radial, or equivalently, it has a monopole-like topology. In this case 
the Poynting flux reads: 



P BH ^--(V/ r B T ,0,0) (30) 



1 

In 

The corresponding electromagnetic luminosity carried out off the black hole reads: 

l bh = J P ™. d S = ±riV v B T B T , (31) 
where ru is the event horizon. 
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Similarly, applying the same procedure to the plasma in the TL, neglecting the effects of vertical advection, and noting that 
Bp = eBj, then the Poynting flux is predominantly radial: 

P TL « = _Ky iB 2^ 0> Q) (32) 

An 

The corresponding luminosity is obtained then by integrating the flux over the radial surface of the TL at r — r tr : 
L TL = j P] h -dS = r tt %V r 4 
Thus, the ratio of the above two luminosities yields: 



(33) 



L BZ 1 V™ Bf n Bj H rii , 

= — (— X- 5 — )(—)(— ) (34) 

How large this ratio is, depends strongly on the location of the transition radius r tt . We now turn to estimate the transition radius 

Noting that at r tr the magnetic energy is in equipartition with the thermal energy of the disk-plasma which, in terms of Eq. A.2, 
means that: 

4 - 05) 
rl p dr 

Furthermore, noting that r BL is the critical radius where B p change topology, and where the corresponding magnetic energy is 
roughly in equipartition with rotational energy, we obtain from Eq. A.2: 

- l T*l (36) 

f r BL P 8r 

Taking into account the spatial variations of B p and p as described in Equations 47 and 50, and solving for X at r BL , we find that: 

C 0.046 if e = 10 1 (— » r tt * 20) 
X(r = r BL ) = e 4/3 ~ I 0.018 if e = 5 x 10~ 2 (— » r tt * 50) (37) 
( 0.002 if e = 10- 2 (^ r tr * 500), 

where r tr is given in r B L-units. Using Eq. 34, and taking into account the spatial variations of Uw and Bj in the TL, and assuming 
that the central BH is rotating at maximum rate, we may obtain the following inequality inequality: 

g < _Ly3/4 ( £H )2= l ( rH )5/4) 
L 1L 4e r tr 4e r tr 

which is larger than 5% for most reasonable values of e(> 10~ 3 ). 

Consequently, the TL appears to provide the bulk of the electromagnetic energy associated with jets, even when the central BH 
is maximally rotating. The contribution of the BHs to the total electromagnetic energy flux is significant, but it is unlikely to be 
dominant. 

Thus, the total power available from the launching layer and from the spinning BH is: 

Lot = L TL + L BZ = L TL (1 + l(r H /r tr ) 5/4 ). (39) 

4e 

If the accretion rate is relatively large and the BH rotates at maximum rate, ttaw, becomes longer than Td yn , and therefore r tr must 
move toward smaller radii. Consequently, the ratio L BZ /L TL increases, but remains below 25% of the total available luminosity 
from the TL for most reasonable values of accretion rates. 

We note that r tl in our present consideration serves as a lower limit. The reason is that the rate of turbulence generation is 
likely to starts decreas ing when the ratio o f the magnetic to internal energies exceeds the saturation value that has been revealed 
by MHD calculations (H awlevetal. 19961) . and is likely to suffer of a complete suppression at r a , where B starts to exceed unity. 

5. Luminosity of truncated disks 

The primary sourc e of heating in standard disks is due to dissipation of turbulence. The corresponding heating function reads: 
( Fran k et al. 19921) : 

1 dQ. , 

D(r) = -W— ) 2 , (40) 

2 or 
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Fig. 8. The distribution of the velocity field in the transition region superposed on equally-spaced isolines of the angular velocity 
Q. The strong-decrease of Q. with radius in the equatorial region relative to its slow-decrease in the TL is obvious. The overlying 
dynamically unstable corona is rotation-free. The orange colour corresponds to large material densities, violet to intermediate and 
blue correspond to extremely tenuous coronal plasma. The 2D distributions of the variable have been obtained using a vertical 
interpolation procedure. 

Please ask A. Hujeirat: 

huj eirat @ mpia-hd . mpg . de 

i ' 

Fig. 9. 20 equally-spaced isolines of the poloidal-component Bp (white contours), superposed on the coloured-distribution of the 
toroidal magnetic field. Blue, green and red colours correspond to low, intermediate and to high TMF-values, respectively. 



where ?7 tur = Sv tur , 2 is surf ace density, and v mr is the standard turbulent viscosity. 

Following the analysis of ( Fra nk et al. 19 92). the r-integrated time-independent equation of angular momentum in one-dimension 
reads: 

3Q , 

ZI/n = 7 7tur — +CI2nr\ (41) 
or 

where C is an integration constant. The classical assumption was that d£l/dr must vanishes at some radius close to the central 
object, and so imposing a physical constrain on C. The idea here was motivated by the two observationally supported assumptions: 
1) Almost all astrophysical objects, that so far have been observed, are found to rotate at sub-Keplerian rates, and 2) v tul is a non- 
vanishing function/process. 

In the recent years, however, two possibilities have emerged: 1) several accretion disk in AGNs appear to truncate at radii much 
larger than their gravitational radii, 2) magnetic fields in accretion disks can be amplified by the Balbus-Hawley instability, 
and potentially they may reach thermal equipartition, beyond which local random motions, turbulence-generation as well as 
turbulence dissipation are suppressed. 

Therefore, we may think of C as some sort of a magnetic effect that come into play whenever the accretion flow approaches a 
certain radius r tt from outside, at which turbulence generation is diminished. In view of Equation 41, we obtain: 

M r tr 3GM M r tr , „ 

77. ur = —[1 - (-) 1/2 ] =» D(r) = -=-3-[l - (-) 1/2 L (42) 
in r iSnr 3 r 

which applies for r > r tt only. The corresponding luminosity reads: 
f" 1 GMM 

L = 2 D{r)2nrdr=- . (43) 

Jt» 2 rtr 

In the case of the nucleus in M87, if the disk surrounding the BH truncates at 100 gravitational radii, then the corresponding 
accretion luminosity is: 

L * ^ = ^ = 2.5 x 10-VXf = 2.5 x lO-VMad * UPergr 1 , (44) 
ZK tr 2<j<J Kg 

where M = 1.6 x 10~ 3 Mem has been used JPi Matteo et a. 20031) . 
Thus, the jet power should of the order of: 

. GMM GMM GMM in44 _, 

L] = * * 10 44 ergs l , (45) 

2r B i 2r BL 2R B l 

which agrees with the observations (see lBicknill & Begelman 19991) . In writing Eq. 45, we have assumed that tbl is located very 
close to the last stable orbit. 
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Consequently, interior to r a the energy goes primarily into powering the jet with magnetic and rotational energies. This power 
dominate the luminosity of the truncated disk, making it difficult to observe the disk directly. 

6. Discussion 

In this paper we have presented a theoretical model for accretion-powered jets in AGNs and microquasars. The model relies 
primarily on electromagnetic extraction of rotational energy from the disk-plasma and the formation of a geometrically thin 
launching layer between the disk and the overlying corona. 

The model is based on self-similar solution for the 3D axi-symmetric radiative two-temperature MHD equations. The very basic 
assumption underlying this model is that standard disks should truncate energetically at a certain radius r tr , below which B p start 
to increase inwards as r~ 2 , and eventually becomes in excess of thermal equipartition. Thus, the matter and the frozen-in PMF can 
be accreted by the BH as far as TAWs are able to efficiently extract rotational energy normal to the disk. Such a flow-configuration 
can be maintained if the PMF change topology at r^L, which must be associated with loss of magnetic flux, so to avoid complete 
termination of accretion. Moreover, unlike standard disks in which the accretion energy is liberated away, or advected inwards 
as entropy (ADAF solutions), the accretion energy here is converted mainly into magnetic energy that power the jet. In this case, 
the plasma in the innermost part of the disk becomes: 1) turbulent-free, 2) rotates sub-Keplerian, and 3) remains relatively cold 
and confined to a geometrically thin disk. On the other hand, the plasma in the TL rotates super-Keplerian, dissipative, two- 
temperature, virial-hot and ion-dominated. 

To first order in e, the profiles of the main variables in the disk are: 



f/ d (r) = [f/2 d+ ^(I_i)]i/2, (46) 

Pdto = P° d T p X-V\ (47) 

Q d (r) = Cl° d X- 5/ \ (48) 

r d (r) = r d °X- 5/8 , (49) 

B p ,dW * B w (r) = B [) A X-\ (50) 



where X = r/r a , T p = [{Ml + lhM^leU^ 12 ^ 2 (see Eq. 17), p° = a£ w Ai U/2 V 15 / 8 f n ' 5 , Q° = r~ 3/2 mdB° = A /^T°. 
Note that the superscript '0' and the subscript 'd' denote the values at the equator and at r a , and which be found using standard 
accretion disk theory. 

In the launching region the variables have the following profiles: 



^(D^d-ilV 4 . (51) 
e a/T 

p w (r) = const. = e —±- U BL , (52) 

r 2 U d -2 V2eA 

n w (r) = Q w T- 5/4 , (53) 

^;w^;, r" 2 , TV) = T w,o Y ~ 5/s < 54 ) 

B p ,wW = B%Y-\ (55) 

fl T ,wto = B%Y- l/ \ (56) 



where Y = r/r BL , A = [(r 7 ^T) - (r 1/4 T) r=ra l Q% = r^ 2 , = 1/e 2 , 7* = 1/e, w = 1/e 2 and B% = l/e 2 . The 
variables here are given in the units listed in Table A.l (see Appendix). 

Furthermore, we would like to point out that: 

1 . Magnetic reconnection in combination with Joule heating is an essential ingredient for our model to work. The magnetic 
diffusivity adopted here, v mag = H w Vj^, is proportional to Vj, and therefore it vanishes at turning points where Bj changes 
its sign. This prescription does not smooth the B T -gradients on the account of Joule-heating the plasma, but it strengthen 
the B T -gradients even more, thereby enhancing the rate of B T -reconnection. The underlying idea here is not to smooth the 
Bp-gradients, but to terminate the propagation of the TAWs to higher latitudes, trapping therefore the angular momentum in 
the TL. In general, a normal magnetic diffusivity cannot hinder the propagation of TAWs into the corona. 
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2. If the axi-symmetric assumption is relaxed, there might be still ways to transport angular momentum from the disk into 
the corona without forming a dissipative TL. Such ways can be explored using a full 3D implicit radiative MHD solver, 
taking into account Bremsstrahlung and Synchrotron coolings, heat conduction, and able to deal with highly stretched mesh 
distributions. Such solvers are not available to date, and therefore beyond the scope of the present paper. On the other hand, 
since the region of interest is located in the vicinity of an axi-symmetric black hole, the axi-symmetric assumption might 
be a valid description for the flow in the TL and in the disk. Furthermore, in the a bsence of hea ting from below, coronae 
around black holes have been found to be dynamically unstable to heat conduction JHuieirat et al. 2002T) . In order to hold a 
BH-corona against dynamical collapse, the MFs threading the corona must be sufficiently strong. These, in turn, force the 
gyrating electrons to cool rapidly through emitting Synchrotron radiation. In the absence of efficient heating mechanisms for 
the plasma in the corona, this cooling may even runaway. 

Consequently, coronae around BHs most likely will diminish, and the robust way to form jets would be through the formation 
of a TL, which might be a natural consequence of accretion flows onto BHs. 

3. We note that if the innermost part of the disk is threaded by an r~ 2 -PMF, then Q ~ r~ 5 / 4 is likely to be the optimal profile for 
the angular velocity. Specifically, 

(a) Q ~ r~ 5 ^ 4 satisfies Eq. 25, so that a mag , and hence M. can be fine-tuned so to produce the radio luminosity and jet-power 
that agrees with observations. 

(b) Multi-dime nsional radiative MHD calculations have shown that r~ 5 ^ 4 is the appropriate profile for the angular velocity in 
the TL JHuieirat et al. 20031) . 

(c) Conservation of angular momentum implies that the flow in the disk must be advective-dominated to be able to supply 
the plasma in the TL with angular momentum, so to maintain the r~ 5 / 4 -profile. 

(d) The Keplerian profile £2 = cr 3 ^ 2 , with c=l yields slow accretion. In this case, the cooling time scales will be faster than 
the hydrodynamical time scale. Thus, the disk cools and the matter settles into a standard disk. However, such a disk is 
unlikely to be stable if the PMFs is in super-equipartition with the thermal energy. 

If c < 1, then r tr must move to much larger radii, and as B p increases inwards as r~ 2 , accretion may terminate by the 
magnetic pressure. 

(e) Adopting the profile £2 ~ r _1 does not fulfil Eq. 25, which requires that the bulk of the generated TMF-energy must be 
advected with the outflow. Moreover, this profile yields unacceptably large and radius-independent toroidal magnetic field 
in the transition layer, and so overestimating the radio luminosity. 

4. The transition radius r tt , which is treated as input parameter, depends strongly on the accretion rate. If the outer disk fails to 
supply its innermost part with a constant accretion rate, then the flow may become time-dependent. In this case, the emerging 
plasmas does not show up as continuous jet, but rather in a blob-like configuration. A high Ai gives rise to ttaw > so 
that r tr must move towards smaller radii. In this case, the total rotational power extracted from the innermost part of the disk 
via TAWs becomes smaller, giving rise thereby to a less energetic and weakly collimated jets. On the other hand, a low VVt 
would allow r tI to move outwards, the total surface from which angular momentum is extracted is larger, and the resulting 
jets are more energetic, have larger magnetic and rotational energies and therefore are more collimated. 

5. We have shown that the BZ-process modifies the total power of the jet. However, we have verified also that the extraction of 
rotational energy from the innermost part of the disk continue to dominate the total power available from the disk-BH system. 

Finally, in a future work, we intend to implement our model to study the jet-disk connection in the microquasar GRS 1915-105 
in the radio galaxy M87. 
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Appendix A: 



Scaling variables 



Mass: 


M- 


= 3 x 10 8 M o 


Accretion rate: 


M 


= l<r 3 Mem 


Distance: 


R = 


Rin = 3R S , where R s = 2GM/C 2 


Temperature: 


f = 


■ 5 x WK 


Velocities: 


V = 


Vs = [yR^T/fr] 1 ' 2 ,^ = 1.23 


Ang. Velocity: 




= V~ Kcp = (GM/R) 1 ' 2 , 6 = V S IV 9 


Magnetic Fields: 


B = 




Density: 


P = 


M/(H d Rl m Vs) = 2.5 x l0- [2 gcm- 3 



Table A.l. The scaling variables used to reformulate the equations in non-dimensional form. 



We use spherical geometry to describe the inflow-outflow configuration around a Schwarzschild black hole. The flow is 
assumed to be 3D axi-symmetric, and the two-temperature description is used to study the energetics of the ion- and electron- 
plasmas. 

In the following we present the set of equations in non-dimensional form using the scaling variables listed in Table A. 1 . 

1 . The continuity equation: 

dp 1 d , 13 

-£ + - T (r 2 P U) + --( COS 0pV) = (A.l) 

ot r l or r cos 8 39 

2. The radial momentum equation: 

dU dU VdU \8P 1 ,1 
ot or r 06 p Or e £ r £ 

1 B e d 113 , 113 , 

+ -^~B X T-(rBi) — (rBl) (A.2) 

p r 00 p r or p r or 

3. The vertical momentum equation: 

dV dV V dV 1 BP 1 V? 
Ot Or r cos 00 rp 00 e 1 r 

_1^_^3 

p r 8r 2rp80 T ' ' 
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4. The angular momentum equation: 

dl 1 d , Id dBj drBg 

+ -j—(rpUi) + --(cosflVf) = Bg— - + B T - 



dt r 2 dr r cos 6 86 86 dr 

cosO 2 8 , dCl 18 , 8Q. 

+ —j-^rirpVtur^r) + a^ COS0 (A.4) 

r z dr dr cos 6 86 86 

5. The internal equation of the ions: 
8&^ 18 18 

-i + --(r 2 P U&t) + --(cosflVfi?) = -(y - l)£fV • V 

8t r 1 dr r cos 6 86 

+ (y - 1)3) - A;_ e + V • K™ nd Vri (A.5) 

6. The internal equation of the electrons: 

-5T + irW) + ^^^ cos0y£ e) - -(r - D£e d v • v 

3f 3r r cos 6 86 

+ (y - 1)0) + A,_ e - A B - A c - A Syn + V • * e cond Vr e (A.6) 

7. The equation of the radial component of the poloidal magnetic field: 

<9f r 2 cos cw 

8. The equation of vertical component of the poloidal magnetic field: 
8Ba 1 8e em r 

-£ H^ =0 ( A - 8 > 

dt rcosd dr 

9. The equation of toroidal magnetic field: 

dBj e d Id 1 dQ 1 dQ. 

— - + - — (rUB T )+-—(VB T ) = -rcos6B r — + -cos6B e — 
dt r dr r d6 e dr e d6 



e z d d Id d 

+ 7^ (Vm ^ rBT) V^ (v ^ BT) (A ' 9) 



where the subscripts "i" and "e" correspond to ions and electrons, p, V = (V r , Vg, V v ) and P are the density, velocity vector, and 
the gas pressure P(= < R Wis p(Tilpi + T e /p. e )), respectively. &f, £f denote the internal energy densities due to ions and electrons, 
& ci = pe,i/(y _ i) f where y = 5 /3 5 ^ = i 23 and p e = 1.14. B = (B r ,B g ,B T ) = (B p ,B T ) is the magnetic field, and e emf is 
the modified electromotive force taking into account the contribution of the magnetic diffusivity. O, A B , A;_ e , Ac, A syn are the 
turbulent dissipation rate, Bremsstrahlung cooling, Coulomb coupling between the ions and electrons, Compton and synchrotron 
coolings, respectively. These processes read: 



$ = v mag |V x B\ 2 /N 



Aj_ e = 5.94 x 10-\iti e ck (Ti - J e) 



3/2 



IN 



A B = 4acK, hsP (T 4 - E)/N, 



A c = 4crn e c( 5 )(T e - T rad )E/N, 

where /c a bs and <x are the absorption and scattering coefficients. N — [(y - l)/y](V 2 V v /R) is a normalization quantity (see Table 
A.l). n e , «i are the electron- and ion-number densities. E is the density of the radiative energy, i.e., the zero-moment of the 
radiative field. The radiative temperature is defined as T ra d = E llA . For describing synchrotron cooling As yn , Eq. 24 is used. 



